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Abstract 

We study invariance and monotonicity properties of Kunita-type sto- 
chastic differential equations in R d with delay. Our first result provides 
sufficient conditions for the invariance of closed subsets of R d . Then we 
present a comparison principle and show that under appropriate condi- 
tions the stochastic delay system considered generates a monotone (order- 
preserving) random dynamical system. Several applications are consid- 
ered. 

Keywords: stochastic delay/functional differential equation, stochastic flow, 
random dynamical system, invariance, monotonicity, random attractor. 

1 Introduction 

In this paper, we study invariance and monotonicity properties of a class of 
stochastic functional differential equations (sfde's) driven by a Kunita-type mar- 
tingale field. Our main results are Theorem 13.41 on deterministic invariant do- 
mains and the comparison principle stated in Theorem 14.21 To prove them we 
represent the sfde as a random fde (see [22] and the references therein). From 
the point of view of deterministic delay systems this random fde has a non- 
standard structure and therefore we cannot apply the results on monotonicity 
available in the deterministic theory. This is why we are forced to develop a new 
method starting from the basic monotonicity ideas. We restrict our attention 
to a class of sfde's which generate a stochastic semi-flow on the state space of 
continuous functions (for an example of an sfde which does not generate such 
a semi-flow, see [19]). For other classes (on L p -type spaces, for instance) we 
can use a variety of approximation procedures to achieve similar results. Our 
choice of continuous functions as a phase space is mainly motivated by the fact 
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that some important results in the theory of monotone systems require a phase 
space with a solid minihedral cone (see, e.g., [Ml fT5] ). 

We note that invariance properties for deterministic functional differential 
equations have been discussed by many authors (see, e.g., [27J[T71[IB] and the 
references therein). We also refer to [55] and to the literature quoted there 
for monotonicity properties of deterministic fde's. Stochastic and random ode's 
were considered in [4]. Similar questions for nonlinear stochastic partial differ- 
ential equations (spde's) were studied in [3J[5] (see also [HI E] and the references 
therein for other applications of monotonicity methods in spde's). 

The paper is organized as follows. 

In Section [5] we introduce basic definitions and hypotheses and describe 
the structure of our stochastic fde model in ([2]) and its random representation 
(see ([5])). The central result in this section is Proposition 12.21 which shows the 
equivalence of the stochastic fde (0) and the random fde ([5]) . 

In Section [3] wc establish our main result concerning the invariance of de- 
terministic domains (see Theorem I3.4[) . The proof involves the random rep- 
resentation established in Proposition 12.21 and also the deterministic approach 
developed in [27j . As an application of Theorem 13.41 we consider an invariant 
regular simplex for stochastic delayed Lotka-Volterra type model. 

In Section 3] wc consider quasi-monotone vector (drift) fields and using the 
same idea as in Section [3] establish in Theorem 14.21 a comparison principle for 
the corresponding sfde's. 

In Section [5] we apply the results of Sections [3J and [4] to construct ran- 
dom dynamical systems (RDS's) defined on invariant regions and generated by 
sfde's from the class considered (see Theorem I5.3|) . These RDS's become order- 
preserving for quasi- monotone drift fields (see Theorem 15 .7|) . In this section 
following pQ (for the monotone case, see also [I]) we recall well-known notions 
of the theory of random dynamical systems including that of a pull-back attrac- 
tor. Theorem 1 5. 71 on the generation of a monotone RDS allows us apply results 
from the theory of monotone RDS's (see, e.g., [2J[U[5] and the literature cited in 
these publications) to describe the qualitative dynamics of the sfde's considered. 
We discuss this issue briefly and provide several examples. 

2 Preliminaries 

Let r > 0, d a positive integer and let C := C([—r, 0], R d ) be the Banach space 
of continuous Revalued functions equipped with the supremum norm || ■ \\c- 
For a continuous Revalued function x defined on some subset of R containing 
the interval [s — r, s] , we define x s G C by 

x s (u) := x(s +u), u € [— r, 0]. 

Let (Q, T, {Tt) t>0 ,P) be a filtered probability space satisfying the usual con- 
ditions. On this probability space we define real- valued random fields M l and 
G% i = 1, 2, . . . , d satisfying the following hypotheses. 
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Hypothesis (M). For each i = 1, 2, . . . , d, M i : [0, oo) x R d x fi ->• R satisfies 

(i) M' is continuous in the first two variables for each we!!. 

(ii) For each x G R d , M*(., x) is a local martingale and M l (0, x, u) = for all 

(iii) There exist S G (0, 1) and predictable processes a IJ : [0, oo) x R 2d x Q — > R 
such that for each i, j G {1, : 



JZ«(t,w):= sup - ' v ' + sup ||D x D„a«(f,s,i,)|| 

^,yeR d I 1 "+" FIJI 1 ~+~ a;,i/eR d 



a ij {t,x,y)\ 

H*l)(i + li 

[jflttfo y) - o'J(t, x, y') - a'J(t, a;', y) + o^t, x\ y') 
X7 t x > tV jiy> \x - x'\ 5 \y - y'\ 5 



is finite, where a y (t,x' ,y') := D x D y a l: > (t,x' ,y') and 



(M I (-,x),M J (-,y)) t = / a l3 (s,x,y,uj)ds a.s., i,j = 1, . . . , d, 
Jo 

where (M l ,M : >) t denotes the corresponding joint quadratic variation (see 
[16] for details). Moreover, we assume that the map t i-> R^(t,uj) is 
locally intcgrable w.r.t. Lebesgue measure for every w S fl and i,j G 
{1, d}. In the definition of D x D y a denotes the matrix formed by 
the corresponding partial derivatives and |.| is an arbitrary norm on the 
space of matrices. 

Hypothesis (G). G = (G 1 , ...,G d ) : [0, oo) x C x f2 -> M d satisfies 

(i) G 1 is jointly continuous in the first two variables for each w£!l. 

(ii) For each u> G f2, bounded set £? in G and T > there exists some L = 
L(T,B,u) < oo such that |G l (i,?7,u;) - G\t, C,cj)| < L||r? - Cllc for all 
< t < T and 77, C G S. 

(iii) For each 77 G G and £ G [0, 00), G(i, 77) is J^-measurablc. 
Below, it will be important to decompose G as 

G(t, v ,uj) =H(t, n, u) + b(t, 77(0), lj), 

where both H and b satisfy (i) , (ii) and (iii) of the previous hypothesis (with 
G replaced by R d with the Euclidean norm for b). In addition, we assume that 
b(t, .) is continuously differentiable for each t and w and there exist 8 > and a 
number c(T, cj) < 00 such that 

sup {su P ||D &( t,x,.)|| + sup \\m^)-m,y,.)\\ , 

0<t<T X £R d x,yeM d ,xjty \ x ~ V\ ' 

(1) 
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In this case we say that Hypothesis (G) holds with decomposition G = H + b. 

For the rest of this section, we assume both hypotheses (M) and (G) and 
fix a particular decomposition G = H + b as above. 

We consider the following Kunita-type delay stochastic differential equation 

dx z (t) = G i (t,x t )dt + M i (dt,x(t)), i = l,2,...,d, t > s, ( ] 

X s = 7], 

where s > and n is a C-valued J^-measurable random variable. 
For the definition of Kunita-type stochastic integrals 

t 

M' l (du,x(u)), 

for adapted and continuous (or more general) processes x, the reader is referred 
to Kunita's monograph [TB] . Readers who are unwilling to learn Kunita integrals 
(even though they are very natural and easy to deal with objects) can think of 
the special case 

m „t 

M l (t,x):=^2 a ik (s,x)dW k {s), (3) 

where W k , k = 1, ...m, are independent Brownian motions and the a lk are 
(deterministic) functions (satisfying appropriate regularity properties). In this 
case @ reads 

m 

dx l (t) = G l (t,xt)dt + J2ci lk {t,x(t))dW k (t), i = l,...,d, t>s, 

k=l 

X s = V: 

and a^(t,x,y) = J2T=i a%k ^i x)a^ k (t, y) is deterministic. 

We aim at a representation of the solution from which one can read off con- 
tinuity properties with respect to the initial condition. Note that even though 
equation © is easily seen to have a unique solution for each fixed s and rj, 
continuity with respect to r\ does not follow since solutions are defined only up 
to a set of measure zero which may depend on 77. To obtain continuity, one has 
to select a particular modification of the solution. We will use a variant of the 
variation- of -constants technique which turns ^ into an equation which does 
not contain any stochastic integral and can therefore be solved for each fixed 
uj G SI. We will see that the modification of the solution which is given by the 
pathwise equation does automatically exhibit continuous dependence upon the 
initial condition. The variation-of-constants technique, which is well-known for 
ode's, has already been applied to sfde's in [50] and [22] . 

For further use we need some properties of the following (non-delay) stochas- 
tic equation 

/ dV/(t) - 6 i (t>(t))dt + M i (dt>(<)), * = l,2,...,d, t>s, . . 
\ ^(s) = x, [ ' 
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where ip = (ip 1 , . . . , ip d ). The following lemma states that equation (Q| generates 
a stochastic flow of diffeomorphisms in M. d . This is a special case of Theorem 
4.6.5 in [H|]. 

Lemma 2.1 We assume that b = (& 1 ,...,6 d ) : [0, oo) x R d x Q ->• M d is a 
vector field satisfying ([1]). T/ien £/iere exists a process '5 : [0, oo) 2 x K' 1 x Q — > M. 
which satisfies the following: 

(i) For each s > and x G M. d , ip(t) = ^> s ^(x,uj), t > s solves equation 

(ii) For each s > 0, x G R d and to G f2, ^ s . s (x,oj) = x. 

(Hi) The maps (s, t, x) <— > \1/ Si t(a;, u>) and (s, t, x) t-» D x ^I' s . t (a;, w) are continuous 
for each uj G fi. Furthermore v f' s .t(.,aj) is a C 1 -diffeomorphism for each 
s,(6l and w G 0. 

flu) For eac/i s,t,u > 0, and w G f2 we Ziawe i/ie semi-flow property 
* S , U {;UJ) = o* Sjt (.,w). 

Note that by (ii) and (iv), we have ^ s , t (.w) = (\& t)S (.,a;)) . 

Lemma 12.11 allows us to construct the following representation for solutions 
to ([2]). In the special case in which the martingale field M is given by a fi- 
nite number of Brownian motions as in ([3]) and 6 = 0, this representation was 
established in Lemma 2.3 in |22j . 

Let * (it, x, u) :=~$ 0tU (x,u). We define the functions £ : [0, oo)xR d xfl -> K d 
and F : [0, oo) x R d x C* x fl -> K by 

a:, a;) := •, w) _1 (a;) = *„,o(a;, w), 

F(u,x,ri,u)) := {D x ^f(u, x, w)} -1 if(u, 77, w) 

and consider the (random) equation 

a;(t, w) = * (t, £(s,r](Q),u) + J F(u,£(u,x(u,w),(j),x u {u),u)du ,ujj (5) 

for t > s with the initial data 

x(t 7 uj) = rj(t — s) for t G [s — r, s], (6) 

where n is a C-valued ^-measurable random variable. We suppress the depen- 
dence of x = (x 1 , . . . ,x d ) on s and r\ for notational simplicity. The following 
proposition shows that equations © and ([5]) (together with ©) are equivalent. 



Proposition 2.2 fix s > 0, a C-valued T s -measurable random variable r\ and 
a stopping time T > s. An adapted M. d -valued process x(t) with continuous paths 
solves equation ([2]) on the interval [s,T(u))]n[s, 00) with initial condition x s = r\ 
if and only if x satisfies ([5]) and ^ on the same interval for almost all lo G fl. 
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Proof. The proof is essentially the same as that of Lemma 2.3 in [52]. Our 
assumptions are slightly different from the ones in [22] but this does not affect 
the arguments in the proof. Therefore, we skip some details. 

First assume that x solves (O and ^ on [s,T(uj)} D [s,oo) for almost all 
w G SI. Then x s = rj almost surely. Equation J5]), together with a slight 
modification of the generalized Ito's formula as stated in [IB], Theorems 3.3.1 
and 3.3.3(i) imply, that a; is a continuous semimartingale and satisfies 



dx(t) =D x -$> ft, [£(s, 7/(0), uj) + F(u, £(u, x(u, uj), uj), x u (uj),uj)du 
x F(t,£(t,x(t),u),x t ,w)dt 

+ ^(dt, [£(s,r/(0),w) + / F(u,£(u,x(u,uj),uj),x u (uj),uj)du ,u) 

J s 

=H(t, x t ) dt + b(t, x(t))dt + M (dt, x(t)). 

Therefore x solves ([2]). 

Conversely, suppose that x solves ([U and define 

C(t,u}):=€(s,r}(0),u})+ / F(u,£(u,x(u),uj),x u ,uj)du. 

J s 

Let 

One can see that x l (t,uj) is a semimartingale with differential 

di'(rj) = H t (t,x t )dt + b l (t,x{t))dt + AP(dt,x(t)), i = l,...,d. 

This (non-retarded) sde has a unique solution x with initial condition x(s) = 
77(0) , so x and x agree on [s — r, T] H [s, 00) almost surely. This proves the 
proposition. □ 

The following proposition provides a well-posedness result concerning prob- 
lem j5j and ©. 

Proposition 2.3 Lei hypotheses (M) and (G) fee satisfied. Then there exists a 
set Slo of full measure such that for all uj £ Qq, s > 0, and r] £ C , the problem ^ 
and ([5J has a unique local solution x(s,T), t, uj) up to an explosion time r(s, 77, uj) . 
For each s > 0, the solution depends continuously upon (t,ij) (up to explosion). 
Further, the following semi-flow property holds: for all < s < t < u, all 77 G C 
and all uj £ Slo, we have 

x(s, 77, u, uj) = x(t, Xt(s, rj, t, uj), u, uj) up to explosion. 

Proof. This is (essentially) Theorem 2.1 in [22]. The only differences are the 
fact that in [22] the authors use the Hilbcrt space M2 instead of C as the state 
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space and that we separate b from G and combine it with the martingale part. 
The proof of local existence, uniqueness and continuity of the problem ([5]) and 
^ is based on a rather standard fixed point argument. □ 

It is natural to ask for sufficient conditions for the explosion time r{s,r],ui) 
to be infinite on a set of full measure which does not depend on s and 77. We 
will say that condition (GE) (for global existence) holds if (G) and (M) hold 
with decomposition G = H + b and there exists a set f2o of full measure such 
that t(s,t],uj) = 00 for all s > 0, all 77 £ C and all u> £ fla- Various sufficient 
conditions for (GE) are formulated in Theorem 3.1 in [22]. They are based on 
spatial estimates on the growth of the flow ^ and its spatial derivative which 
were established in [21] and [12]. We quote them here: 

Proposition 2.4 Let (G) and (M) hold with decomposition G = H + b. Each 
of the following conditions is sufficient for (GE): 

(i) For each T > and u E Q there exist c = c(T, u;) and 7 = 7(T, w) £ [0, 1) 
such that 

\H(t,r,,u)\<c(l+\\r,\\l) (7) 
for all Q <t <T , ^ £ C and to £ f2. 

(ii) For each T > there exists /3 £ (0,r) such that H(u, ?/,w) = H(u,fj,uj) 
holds for all to E £1 whenever < u < T and r)\\_ rt _p] = ri\[- r -/}]• 

(Hi) For all lo e fl and T G (0, 00) we have that 

sup ||(D a; ?/i(u,x,a;)) _1 || < 00 

0<u<7>GR d 

and there exists c = c{T,uj) such that ([7]) holds with 7 = 1. 

It is a bit annoying that (i) excludes the case of H satisfying a global Lips- 
chitz condition. It seems to be open whether (GE) holds in that case. 

3 Deterministic invariant regions 

In this section we assume that Hypotheses (M) and (G) with decomposition 
G = H + b and condition ([7]) are in force and consider a general problem of the 
form (O, ([6]). We provide sufficient conditions that, given a non-empty closed 
(deterministic) subset D in R d , a solution with values in this set for t £ [to — r, to] 
will have values in D for all t > to- The key idea is to decompose the solution 
semi-flow in such a way that ^ alone leaves D invariant and that the remaining 
drift does not change this property. 
Below we use the notation 

C = {r} £ C : r)(s) £ D for every s £ [-r, 0]} . (8) 
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We need some additional hypotheses (which are inspired by similar hypotheses 
for deterministic fde's in [27]). 

Hypothesis (G £ ). There exists a family {G e } of random fields satisfying (G) 
with decomposition G e = H £ + b for every e £ (0, £o] such that 

(i) H e satisfies condition ([7]); 

(ii) lim e _,.o H e (t, tj, oS) = H(t, rj, w) for every (t, rj) £ [0, oo) x Cd and uj £ CI; 

(iii) given (t, if) £ [0, oo) xCd and e £ (0, £q] there exists an a = a(e, t, t),uj) > 
such that if < h < a and u £ M. d is such that \u\ < a, then 

r)(0) + hH e (t,r),u)) + hu £ B; 

(iv) if y e (t, rj) solves the problem (JSJl,© with G £ instead of G, then for every 
uj £ CI, t > 0, and rj £ C we have lim e ^o y £ (^ ^ 5 w ) = %(t,r),u)), where 
x(t,rj,uj) solves (JS]),©. 

Note that condition (ii) implies that D is the closure of its interior. Further 
note that condition (iv) above is not implied by the other conditions - not even 
in the case of a deterministic ode, see [27] . 

Hypothesis (Mb). The problem (|4]) generates a stochastic flow $ t)S (-,u>) of 
diffcomorphisms of K d such that 

* a ,t(B,a;) =B, t > s, uj e tt. (9) 

Remark 3.1 For flows which are driven by a finite number of Brownian mo- 
tions, explicit criteria for the validity of this hypothesis are well-known (wc 
will state some of them below) . We have not found corresponding criteria for 
Kunita-typc equations in the mathematical literature. In fact, such criteria fol- 
low easily in case B is compact: for Kunita-type sde's, the one-point motion 
(i.e. the solution for a single starting point x £ M d ) can be described by an 
equivalent sde which is driven by a finite number of Brownian motions (which 
depend on the point x). Assuming that for each x £ B the solution stays in B 
forever with probability one (for which one can check the known criteria), then 
the same holds true for a countable dense set of initial conditions in B. The 
fact that B is compact and the flow is continuous shows that there exists a set 
f2o of full measure such that © holds for all ui £ Q. Our claim follows since we 
are free to modify $ona set of measure zero. 

Remark 3.2 If B is a closed convex subset of M. d with nonempty interior then 
Hypothesis (G e ) follows from the Nagumo type relation 

lim + /i" 1 dist(ry(0) + hH(t, ?], w), B) = (10) 

for any t £ M, w £ CI and 77 £ C such that 77(5) £ B for s £ [— r, 0]. In this case 
we can take 

H e (t, r), u) = H(t, 77, w) - £(77(0) - e), 
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where e is an element from hit D. If D = R 1 ^ relation (|10p is equivalent to the 
requirement 

{n > 0, J7*(0) = 0} => H l {t,r),uj) > (11) 
for every t, u> E and i = 1, . . . , d. For the proofs we refer to [27] . 

In the following remark we discuss conditions and examples when Hypothesis 
(Mb) is valid. 

Remark 3.3 Assume that D is a closed set in R d such that D has an outer 
normal at every point of its boundary. We recall that a unit vector v is said to 
be an outer normal to D at the point xq € <9B, if there exists a ball B(x\) with 
center at x\ such that B{x\) flD = {xq} and v = X ■ (x± — xq) for some positive 
A. 

Let W\, Wi be independent standard Wiener processes. We consider ((4]) 
with 

i 

M\dt,m) = x>j(<K*)w-(t), (12) 

where the coefficients have bounded derivatives up to second order. The problem 
in @ can be written as a Stratonovich sde: 

i 

dip l (t) = 6 i (i,V(i))df + £ro*-(^(f)) °dWj(t), i = l,...,d, 

-0 l (s) = X, 



where "o" denotes Stratonovich integration and 



l d 



dm l Ax) 



^ ' dx k ' 

3=Xk=l 

It follows from Wong-Zakai type arguments that D is forward invariant under 
tp if for any x € 9D we have 

d d 

^6 4 (i,a;)^<0 and ^m*(x)4=0, j = l,...,l, (13) 

i = l i=l 

for every outer normal v x = (i/^, . . . , i/^) to D at x. We refer to [4j Chap. 2, 
Corollaries 2.5.1 and 2.5.2 ] for details. Further, (Mj) holds if D is both forward 
and backward invariant under ip. For this to hold it is sufficient to assume that 
the first inequality in ([TB")) is an equality for each x € 3D. We note that in the 
case D = and W(t, x) = the first condition in (| 13[) follows from the second 
one which can be written in the form 

m)(x) = for all x = (xi, . . . , Xi-i, 0, Xi+i . . . , Xd), i = 1, . • . , d, j = 1, . . . , I. 



As an example we point out the case when D = {(x\; x-i) : x\ + x\ < 1} Cl 2 
and problem ((4]) has the form 

dxt = m l (xi,x 2 )dW(t), i = 1, 2. 

In this case b l (t,x) = and relations (|T3|) holds if ra (ai) = m 2 (x) = for all 
\x\ = 1. For instance, we can take m 1 (xi, X2) = — c(|x|)x2 and m 2 (xi, X2) = 
c(|x|).ti, where c(r) = for r = 1. 

Our main result in this section is the following theorem. 

Theorem 3.4 Assume that Hypotheses (G), (M) ; (Mb) and (G £ ) hold. Then 
ID) is a forward invariant set for problem @ m i/ie sense that for any uj G f2, 
s > and n £ C such that rj{u) G D /or u G [— r, 0], £/ie (unique) solution x of 
(|5|) and ((HJ) satisfies x{t,r],ui) G D /or aZZ t > s. 

Proof. Assume that D is not forward invariant. Then there exist w G f2, s > 0, 
77 G Cd and t* > s such that x(£*) D. Let y e (i) be a solution to the auxiliary 
problem in (G e )(iv). It follows from assumption (iv) in (G e ) that there exist 
e > and to G [s, £*) such that 

y £ (t) G D, t G [s - r, t„] and j/ e (t„ + fy) g D, 

where {/ij} is a sequence of positive numbers such that lim :( _s. 0O hj — 0. The 
solution y e (t) can be represented in the form 

y e {t,w)=i>{t,C, e {t,u 1 ),u>), t>t , 

where 

( £ (t,uj) = Z(t ,y £ (t ),uj) + f F £ (u,^u,y £ {u),uj),y £ u ,uj)du. 

Jto 

Here £(t,x,uj) = ip(t, ■,uj)~ 1 x, ip(t,-,ui) — Wo,* (•, <^), where ^o,t(-,^) is the 
diffeomorphism given by (UJ and 

F £ (u,x,ri,uj) = {D x i/)(u,x,cj)}~ H £ (u,n,uj). 
Since (u, x, 77) 1 — y F £ (u, x, 77, ui) is continuous, we have that 

rto+h 

/ F £ (u^(u,y £ (u),uj),y £ ,,uj)du = hF £ (t ,£(t ,y £ (t ),uj),y £ o ,uj) + o(h). 
Jt 

Thus 

C £ (t Q + h,w) = £(t ,y £ (t ),Lj) + hF £ (t ,Z(t ,y £ (t ),Lo),y £ o ,uj) +o(h). 
We have that ip(to, ^ _1 (io, x, uj), uj) = x. Therefore by the chain rule 
DzV^o, ?/'" 1 (to,x,a;),a;)D x '0" 1 (io,a;,iu) = Id. 
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Thus 

{D z V(io,^ _1 (io,a;,^),w)} = D x V> _1 (io,£,w). 

Consequently, 

F E (t , £(i , y e (to) , w) , vt , w) = D, i>~ 1 (t , 2/ £ (to) , wjfle (to , y! , w) . 
It is also clear that 

IT 1 (*o , V s (h ) + hH E (t , y! o , w) , «) - V" 1 (to , V s (to ) , w ) 

= hV x ^-\to,y s (to),w)H e (to,y s t ,u) + o(h). 

Thus 

C e (t + h, u) = V-" 1 ^, y e (t ) + hH E (t , y e tQ , uj), u) + o(h). 
This implies that 

^(to,C(to + h,u),u) = y E (to) + hH £ (t ,yl ,uj) + o(h). (14) 

Hypothesis (G E )(iii) implies that the right-hand side of (TT4")) is in B for all 
sufficiently small h > 0. By Hypothesis (Mb) wc therefore have ( e (to+h, u)eD 
and hence y e (io + h, uj) = ip(to + h, ( e (to + h, uj), uj) E B for all sufficiently small 
h > contradicting our assumption that y e (to + hj,ui) £ B for all j. This 
contradiction proves the theorem. □ 

In the following assertion we show that, similarly to the deterministic sit- 
uation (see [2S] and the references therein), in some cases the Nagumo type 
condition in (|10p provides us necessary and sufficient conditions for invariance. 

Corollary 3.5 Let B be a closed convex subset of K d with nonempty interior 
and let Hypotheses (G) , (M) and (Mb) be in force. Then B is a forward 
invariant set if and only if (|10[) holds. 

Proof. If (fit))) holds, then we can apply Remark 13.21 to conclude that B is 
forward invariant. 

Let B be forward invariant. It is clear that fT4")) holds for e = and to = s, 
i.e. we have 

ip(s,C(s + h,u),u) = r)(0) + hH(s,r),u) + o(h) (15) 

for any rj E Co- Since B is invariant, we have tp(s + h, C(s + h, uj), uj) = x(s + 
h, r],uj) E B for all h > 0. (Mb) implies that ip(s, £(s + h, oj),oj) lies in B for all 
h > 0. Therefore (JTSJ) implies (HUJ). □ 
In the case B = Rl Corollary 13.51 implies the following assertion. 

Corollary 3.6 Let (G) and (M) fee in force with M 2 of the form (| 12[) . Assume 
that 

b l (t,x) = 0, m}(a;) = V a; = (xi, . . . , Xi_i,0, . . . , Xd), 

where i = 1, . . . , d, j = 1, . . . , I. Then is forward invariant set if and only if 
([XT]) ZioZds. 
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Proof. It follows from Corollary 13.51 see also Remarks 13.31 and 13.21 □ 

More complicated example of an invariant set D is discussed in the following 
remark. 

Remark 3.7 Assume that D is a set of the form 

D = {x 6 K d : {a q ,x) < 7 „ q = 1, . . . , Q} , 

where a q = (a q , . . . , a d ) <E R d , j q & q = 1, . . . ,Q. Then ([TO) holds if and only 
if for every q = 1, . . . , Q we have the relation 

d 
i=l 

whenever rj £ Co and (a q ,r)(0)} = j q . 

In the case when b % = and M l has the form (jT2"j) . it follows from Remark [3. 3 1 
that Condition (Md) holds if for every g = l,...,Qwe have the relations 

1 d d dmUx) d 

j — l k—l i—1 2—1 

whenever {a q ,x) = 7 9 . For instance, this condition is true if 

Q 

m){x) =a t j (x 1 ,...,x d ) ■ Y[h q ((a qi x) -7,), 

q=l 

where crj(xi, . . . ,Xd) are arbitrary and h q (s) is such that h q (0) = 0. 
Now we provide some examples. 
Example 3.8 For the system 

1 

dx\t) = x l {t)f\x t ) dt + x l {t) ^(^(t), ...,x d (t))dW j (t), i = l,...,d, 

i=i 

the set IR+ is a forward invariant set. Here /' and cr* are such that condition 
(GE) holds. This conclusion follows from Corollary 13.61 

Example 3.9 In the previous example, the noise can also be replaced by a 
more general Kunita-type noise. As a particular example, let N be space-time 
white noise on M. d x [0,oo), let h : M. d — >• [0,oo) be C°° with compact support 
and define 

M\dt, x) := 4>{x l ) [ h(x - z)N(dz, dt) 7 
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where <j> G C°° is bounded and all its derivatives are bounded. Assume that 
0(0) = 0. Then 

{M%,x),MS(.,y)) t = <Kx i )<l>(y1)t f h(x - z)h(y - z)dz. 

Note that M satisfies hypothesis (M) due to our assumptions on and h. To 
facilitate things, we assume that f l (r)) = g l (f](—r)) where g % is bounded and 
Lipschitz. Then the set R^ is a forward invariant set for the flow generated by 

dx*(t) = ^(t^g^xit - r)) dt + M\dt, x(t)), i = 1, d. 

To see this, first note that condition (GE) holds by Proposition 12. 4( i). Then 
we argue as in Remark 1 3. II for each given starting point x G R+, the solution 
starting in x will remain in R^ forever almost surely since it can be written as 
a solution to an sfde with finitely many driving Brownian motions, so the same 
property holds for all starting points in R^J. with rational coordinates. Since we 
know that a local flow exists, no trajectory of the flow can leave Rl through 
one of the hypcrplancs bordering Rl. Since also (GE) holds, no trajectory of 
the flow can escape to infinity at finite time either, so the set R+ is invariant. 

Note that in this set-up the driving noise M is independent at locations x 
and y with distance larger than the diameter of the support of h which is a 
reasonable assumption in many models and which cannot be achieved with a 
finite number of driving Wiener processes. 

Example 3.10 (Lotka-Volterra type model) Consider the system 

dx\t) = -aiX l {t){l-{b,x{t-r)))dt + aiX l {t){l^{b,x{t)))dWi{t), i = l,...,d, 

(16) 

where b G K+, on > 0, cr t G R. The set D = {x G R+ : (b,x) < l} is forward 
invariant. Since D is a bounded set in R d . we can modify the nonlinear terms 
outside some vicinity of ID in order satisfy the requirement in ||7J) This allows us 
to apply Propositions 12.31 and I2.4f i) and obtain well-poscdness of the problem 
in (|16[) . The statement on the invariance follows from Theorem 13.41 via the 
observation made in Remark 13.71 

4 Comparison theorem for sfde's 

Our next result is a comparison principle for functional differential equations 
perturbed by Kunita type noise of the form ^ with the local martingales M l 
not only satisfying Hypothesis (M) but also 

M l (t, x,(j) = M ! (t, Xi,uj) for all i=l,...,d, x = (x\, . . . , Xd), 

i.e., M % depends on t, u> and on the i-th component of spatial variable x only. 
As can be seen from [4] , this structural requirement is needed for a comparison 
principle even in the non-delay case. 
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Thus, instead of @, we consider the following Kunita-type retarded stochas- 
tic differential equation 

dx l (t) = G l {t,x t )dt + AP{dt,x l (t)), i = l,2,...,d, t>s, , 
%s = V, 

where 77 is a C-valued J^-mcasurablc random variable and the drift terms G l 
satisfies (G). We fix a decomposition G = H + b as in the previous section 
and assume that b l depends on x l only. In this case instead of ((4]) we have the 
diagonal system of scalar non-delay equations 

<ty*(f) = V(t,^{t))dt + NP{dt,^ l {t)), i = l,2,...,d, t>s, 

It follows from Lemma 12.11 that equation (fT5|) generates a stochastic flow of 
diffcomorphisms x h-> ip z s t (x, ui) in R for each i = 1,2, ... ,d. Moreover, 

V s ,t(x,u>) = (^ )t Oi,Cj), . . . ,ij) % Stt {x d ,u))) 

satisfies all statements of Lemma \2. II Below we often write ip l (t,x,ui) instead 
of ipQ t (x,ui). Observe that due to the diffeomorphic property the flow ip l is 
automatically strongly monotone in the sense that for x, y € R we have 

x < y implies ip l s t (x,uj) < ip l st (y,oj) for each s,t > 0, u € 0. 

Indeed, if the implication above is not true, then there exist x < y, s < t, and u 
such that ipl t (x,uj) = ip\ t {y,ui). Since ip l s t (-,uj) is invertible, this implies x = y 
and thus provides a contradiction. 

Applying Proposition 12.21 we can specify representations (|5|) and © for our 
case of diagonal M 1 . Namely, if we define 

C(u,x\uj) := ^(tv.wrV) 
F i (u,x i ,T),u) := {0^(11, x^u)}' 1 77, w), 

then ([5]) can be written in the form 

x i (t,u)=^ i (t, ?(a,ri i (0),u) + J F\u,C{u,x\u,oj),oj),x u {uj),uj)du ,u 

(19) 

for all t > s. 

Let C+ be the standard cone in C. This cone defines a partial order relation 

via 

77 > 77, iff r\ - 77* € C+, (20) 
i.e., iff 77* (s) > 77* (s) for all s g [— r, 0] and i = 1, . . . , d, where 

ri= (77 1 ,...,7 ? d ) and r)* = (vl,---,vt) 
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are elements from C = C([— r, 0], M. d ). We write rj > 77* iff 77 > 77* and 77 7^ 77* 
and use the notation 77 >> 77* if 

?f (s) > 7/* (s) for all s G [— r, 0] and i = 1, . . . , d. 

We also consider another sfde 

f cLc^i) = G i (<,x t )d7j + M i (dt,a; l (t)), i = 1,2, . . . ,d, t > s, ( . 
\ x s = 77* e C, ^ ^ 

with the same M and 6. We assume that the random field G = {G 1 } satisfies 
Hypothesis (G) with decomposition G = H + b. Let B(w) C K d be a closed set 
with nonempty interior. 

Definition 4.1 Let D(w) C K d be a closed set with nonempty interior. For 
each uj, let [a(u),b(uj)] a random interval in R + . A random vector field G = 
(G 1 ,...,G d ) : [0, 00) x Cx Q ->• R d is said to be quasimonotone on [a(w), x 
D C M^ 1 iff for any 7/ = (^(s), r? d (s)) and 77* = (?yj( s )> • • • ,?7*(s)) from C D , 
where Cd is defined by ([5]), we have the following implication 

{v > V*, V l (0) = Vl(0)} => G\t,n,u) > Gifaw) 
for every t G [a(w), b(oj)], uj G £1 and i = 1, . . . , d. 

We note for future use that quasimonotonicity is invariant with respect to a 
decomposition G = H + b with b % depending on x l only in the sense that G is 
quasimonotone if and only if H is quasimonotone. 

Theorem 4.2 (Comparison Principle) Assume that M satisfies the condi- 
tions above and that the random vector fields G = (G 1 , . . . , G d ) and G = 
(G 1 , . . . , G d ) satisfy Hypothesis (G). Let x(t) := x{t, r),u>) be a solution to (|17|) 
and y(t) = y(t,n,uj) be a solution to ([2"Tj) which possess the property 

x(t), y(t)eB, for t G [s, s + T(uj)} 

for some convex closed set D C R d with nonempty interior, where T(ui) > for 
all uj G Q. Assume that (Mi) holds and the random field G is quasimonotone 
on [s, s + T(lu)] x ES. Then the following assertions hold: 

!• If V ^ V* an d 

G(t,t,u)<G(t,£,w) for all £ G Co, t G [s, s + T(oj)}, uj € fi, (22) 
then 

x(t; 77, uj) < y(t, 77*, uj) for all t G [s, s + T(uj)], uj G fl. (23) 

>G(t,£,w) for all £ G C D , f G [a, s + T(o>)], w G fi, (24) 
ir(t; 77, w) > j/(t, 77,, w) /or a/Z i 6 [s, s + T(w)], a; G fi. (25) 
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Proof. We prove the first part only (the proof of the reversed inequalities is 
similar) . 

We start with the case rj << 77* and G « G, i.e., we assume that 

V l {s) < Vl(s) for all s G [-r, 0] and i=l,...,d, (26) 

and 

G'(t^,uj) < &(t,^u) for all £ G C®, * G [s, s + T(u)}, u G Q, (27) 

where i = 1, . . . , d. The same is true for H and H Let us prove that 

x i (t;r],uj) < y\t,r]*,u)) for all i G [s, s + T(u>)], w G 0, i = 1, . . . , d. (28) 

Since x(t) and y(i) are continuous for all u G fi relation (|28[) is valid for some 
interval [s,s + t(ui)], where < t(uj) < T(lo). If (|2"5)l does not hold for a// t 
from [s, s + T(w)], then for some ui there exist t' G (0, T(w)) and i G {1, . . . , d} 
such that 

x l {t') = y l (t') and ar'(i) < y j {t) for all t G [s,s + i') 5 J = 1, ■ • ■ ,d. 
Using representation (|19[) and strict monotonicity of ip 3 we obtain that 

Ch(t') = Ch(t') and CM*) < C&(*) for a11 * e [«> s + 0. J = 1, • • • > (29) 
where 

C G (*) := ? (s,rf (0), w) + j F 3 G (u, ? (u, x 3 (u),u),x u ,u)du, 

C&(t) :=?(a,ri(0),w)+J F 3 Q {u,£, 3 \u,y 3 \u),w),y u ,u)du 

with the following notation: 

£ J (u, x J , w) := ip 3 (u,-,oj)~(x 3 ), 
F G (u,x 3 ,rj,u) := {D^ip 3 (u,x 3 1 H J (u,t] 7 uj) 7 

Fq(u,x? ,t/,w) := (D !E j^ 3 (u,x , ',a;)} 1 H 3 (u,t],uj). 

Since the functions F G (u, x 3 , rj,oj) and Fjk(u, x 3 , rj, w) are continuous for every 

cj G f2, the pre 
the equations 

dt' 



G f2, the processes CgW an d CgW are continuously differentiable and satisfy 



-e G (t) = F 3 G (t,t?{t,x 3 \t),u),x t ,u) (30) 



and 

d 

dt 



^(*)=-FS(t,^(*,y , '(*),w),» t ,w). (31) 
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It follows from that 

Chit') - Ca(t) > Chit') ~ Qf(*) for all s < t < s + t'. 
This implies that 

^(0 > ^(f). (32) 
However, since = y l (t'), from (|3Tj) we have that 

^CM*') = {Dx^'Ka: i (0^)} _1 lf i (t ! ^,a;). 

Therefore (|2"T)) written for if and H, quasimonotonicity of H and ([50)) imply 
that 

> {D x ^(« > x i (* , ),w)}" 1 iP(t / ,Xf ,w) = ^CgCO- 



This relation contradicts to ([52]). Thus {55} and (gTJ) imply 

To prove for the general case we first apply the result above to the 
corresponding equations with 

and 

G e (t,ri,w) = G{t,r J ,w)+e{e 3 -r](0)), 
where ei, e2 S intD and ei << e2. It is clear that ([22]) implies that 

G*(t,f,w) <G*(t,£,w) for all £ e G D , t E [s, s + T(w)], cj e fi, * = 1, . . . , d. 



Thus, by limit transition we obtain (|23j) in the case when r\ < < 77* . Using this 
fact it is easy to prove ([23]) for every 77 < 77*. □ 

Remark 4.3 If the drift term G is quasimonotone on R+, then we have that 
G(t, 77, lj) > G(i, 0, w) for every 77 6 G K d . Therefore applying the comparison 

principle in (|2"4"1) and (f23|) with G = we can conclude that is a forward 
invariant set with respect to sfdc (fT7|) when G(t, 0, w) > and M (t, 0, w) = 0. 

Example 4.4 (Lotka-Volterra type model) Consider the system 



dx l {t) = OiX^t) I 1 - fcx^t) - y ^Cij I x j (t + T)dfJ,ij(r) \ dt 



(33) 



+<7iX l {t)(Ri - x l (t))dW u t> s, i = 1, 
77 g G. 
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Here at, Pi and Ri are positive numbers, Cij > 0, &i G K. We assume that Hij(r) 
are left continuous nondecreasing functions on [— r, 0] of bounded variation such 
that 

A%(0) - Vij{-r) =1, i,j = 

It is easy to see from Theorem l3.4l (see also Corollary 13. 61 and Remark l3.7p that 
D = riiLiPj-Ri] is a forward invariant set for sfde (f33|) provided i?i > /J^ 1 for 
every i = 1 , . . . , d. 

We note that the global well-poscdncss of (f3U|) follows from Propositions 12. 31 
and I2.4f ii) because we can modify the corresponding drift term outside D to 
satisfy Q. 

It is also clear that the functions 



&(r,) := a i? /(0) [ 1 - /W(0) - °« f 



rf(T)dmj(r) , i,j = l,...,d, 



satisfy the inequality 



i? f(0) ( 1 - f3rf(0) - J2 °^Rj \ < &{r,) < <W(0) (1 - ptfiO)) 

for every -q G Cm, where Cd is given by © with D = JliLi [0, Ri]- Since the 
functions 



Gjfa) := ^(0) 



and Gjfa) :=cvf(0) (l-ftV(O)) 



are quasimonotone, Theorem 14.21 implies that for any initial data r\ G Cd a 
solution 

x(t,r), cj) = (x 1 ^,^,^), . . . ,x d {t,-q,u)) 
to problem (|33[) satisfies the inequality 

u l (t,r],u}) <x l (t,r},u) <v l (t,r),w), i,j = l,...,d, (34) 

where u(t, rj, uj) = rj, u>), . . . , u d (t, i], us)) solves the problem 



du'(t) = a lU l {t) ( 1 - piU^t) - Y c v R j ) d * 

3=1 



I «*(0) 



+cr i u i (t)(R i - u i (t))dW l , t> s, i = l,...,d, 
min re[ _ r , ] V l (s), i = 1, 



..,d, 



(35) 
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and v(t, rj, w) = r), lo), . . . , v d (t, rj, w)) solves the problem 

( dv l (t) = oiiV^t) (1 -PiV^t)) dt + cTiv'itXRi -u i (t))dW i , t > s, 
[ w 4 (0) = max Te[ _ r . 0] rj l {s), i = l,...,d. 

(36) 

We emphasize that problems ([33)) and ([33)) are direct sums of one-dimensional 
ordinary stochastic differential equations. Long time dynamics of these ID 
systems is described with details (see, e.g., @] and the references therein). Thus 
we can use the relations in (f34|) to "localize" dynamics of the original sfde ([33)) . 

5 Order-preserving RDS generated by sfde's 

In this section we consider some other applications of Theorems 13.41 and 14.21 
from point view of theory of random dynamical systems (RDS). 

5.1 Generation of RDS in an invariant region 

Following the monograph of Arnold [T], we introduce the notion of a random 
dynamical system. 

Definition 5.1 Let X be a topological space. A random dynamical system 
(RDS) with time K+ and state space X is a pair <fi) consisting of the following 
two objects: 

1. A metric dynamical system (MDS) d = (0, T, P, {i9(i), t £ R}), i.e., a 
probability space (f2, J-, P) with a family of measure preserving transfor- 
mations {$(£) : H *-> fi, t £ R} such that 

(a) tf(0) = id, o tf(s) = tf(t + s) for all t, s £ R; 

(b) the map (i, w) h-> $(t)uj is measurable and #(i)P = P for all 

2. A (perfect) cocycle </> over $ of continuous mappings of X with one-sided 
time R+, i.e. a measurable mapping 

: R + x x X h-> X, (t,iM,x) <-> <fi(t,uj)x 

such that (a) the mapping cf>(-,uj) : x n> <j){t,w)x is continuous for all 
t > and wefi; (b) it satisfies the cocycle property: 

</>(0, w) = id, <f>{t + s,ui) = <j>(t, ■d(s)uj) o cj)(s, uj) 

for all t, s > and u> £ 0. 

Definition 5.2 Let $ be an MDS, J" the P-completion of T and F = {J"t,t G 
R} a family of sub-cr-algebras of T such that 

1. T s C Ji, s < i; 
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2. T s = f\>o Fs+h, s£R, i.e. the filtration F is right-continuous; 

3. J- s contains all sets in T of P-measure 0, s € 1; 

4. i?(s) is (Tt+s, J r t)-measurable for all s,( £l. 

Then (i), F) is called a, filtered metric dynamical system (FMDS). If - in addition 
- (■#, <fi) is an RDS such that (f>(t, -)x is ( Ft, £>(X))-measurable for every t>0,i£ 
X, then ($,F, <fi) is called a filtered random dynamical system (FRDS). 

We recall that an X-valued stochastic process Y(t),t G T C R is called adapted 
or nonanticipating with respect to the filtration F if is (J-" t ,£>(X))- mea- 
surable for every t G T. Therefore (0,F,0) is an FRDS iff is an RDS, 
($, F) is an FMDS and </>(•, -)x is adapted to F for every x G X. 

Theorem 5.3 Assume that Hypotheses (Mb) and (G e ) are in force and (GE) 
(see Proposition \2.J$ holds. If the drift term G(t,r],uj) = G(q) and the the local 
characteristic a of M are deterministic and autonomous, then problem ([5j • © 
(and hence @ ) generates a FRDS ($, <p) in Co, where Co is defined by (|SJ| (the 
case © = R d is not excluded). The corresponding cocycle tp has the form 



[<p(t,w)r,](T) 



x(t + T,T),u), t + T>0, 

V(t + r), t + r<0, 



for every r G [— r, 0], where x(t,rj,w) is a solution to problem ^ with s = 0. 
Moreover <p(t,w) is compact mapping in Co, i-e. for any bounded set A from 
Co, the set <p(t,u)A is relatively compact in Co for every t > 0. 

Proof. It follows from Theorem 13.41 and from the representation in ([S]),© of 
solutions to (0) . We also use Propositions 12.21 12.31 and 12.41 □ 

To describe long-time dynamics of an RDS we need a notion of a random 
set (see, e.g., [T] and the references therein). 

Definition 5.4 A mapping cj h-> D(uj) from H into the collection of all subsets 
of a separable Banach space V is said to be random closed set, iff D(uj) is a 
closed set for any well and w i— > disty(x, D{uj)) is measurable for any x G V . 
The random closed set D(uS) is said to be compact, if D(ui) is compact for each 
uj. The random closed set D(lo) is said to be tempered if 

D(oj) C {x G V : \\x\\ v < r(uj)}, lo G O, 

where the random variable r(ui) possesses the property sup tgK {r(??(t)a;)e _7 '*' } < 
oo for any 7 > 0. 

We also need the following concept of a random attractor of an RDS (see 
[HI US] and also [J 2] and the references therein). Below we denote by X any 
subset af separable Banach space V equipped with the induced topology. 

Let V be a family of random closed sets in X which is closed with respect 
to inclusions (i.e. if D\ G T> and a random closed set {£>2(w)} possesses the 
property D%(u)) C Di(u>) for all to G f2, then D2 G £>). Sometimes the collection 
V is called a universe of sets (see [T]). 
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Definition 5.5 Suppose that ("#, </?) is an RDS in X. Let V be a universe. A 
random closed set from T> is said to be a random pull-back attractor of 

the RDS (/?) in V if 7^ X for every w £ f2 and the following properties 
hold: 

(i) A is an invariant set, i.e. (p(t, lj) A(u>) = A($(t)(jj) for f > and w£(l; 

(ii) A is attracting in 2?, i.e. for all D 6 V 

lim cMipft, t?(-*)w)D(i9(-i)w) I A(lu)} = 0, w £ ft , (37) 

t— »+oo 

where dx{A\B} = sup x£A distx{x, B). 
If instead of (|37l) we have that 

lim P{w : dx{<^(t,w)-D(a;)| J 4(^(f)a;)}>5} = 

i— 5- + 00 

for any 5 > 0, then is said to be a random weak attractor of the RDS 

Some authors (e.g. [HI H]) require a random attractor to be compact (and do 
not insist that it is different from the whole space). This distinction will not be 
important in what follows. The notion of a weak random attractor was intro- 
duced in |23j . For the relation between weak, pull-back and forward attractors 
we refer to [25]. 

Remark 5.6 If D is bounded in M. d , then it is easy to see that the RDS (1?, ip) 
generated by ^ in Co has a random compact pull-back attractor in the universe 
V all bounded sets. Since by Theorem 15.31 the RDS ($, ip) is compact, this 
follows from Theorem 1.8.1 [1], for instance. In the case of unbounded sets D 
(e.g., D = Ml) we need some conditions which guarantee dissipativity of the 
corresponding RDS. These conditions can be obtained in the same way as for 
the non-delay case (see, for instance, [H Theorem 6.5.1]). 

5.2 Monotone RDS 

Let as above C = C([— r, 0],K d ) and C+ be the standard cone in C of nonneg- 
ative elements: 

C+ = {r ] = (v 1 (s),...,r ] d (s))eC V i (s)>0 Vse[-r,0], i = l,...,d}. 

This cone is a normal solid minihedral cone. This fact is important for further 
application of the theory of monotone RDS. We refer to [14] and [15] for more 
details concerning cones and partially ordered spaces. 

Let D be a convex closed set in R d with nonempty interior. In the space 
Co given by ((SJ) we we define a partial order relation via (|20p . i.e., r\ > 77* iff 

t] - e C+. 
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Theorem 5.7 Assume that the hypotheses of Theorem \ 5.S\ hold. Let D be a 
convex closed set in M. d with nonempty interior and tp) be the FRDS generated 
by problem (|17[) in Co defined by ([5]). If the random field G is quasimonotone 
in D then 93) is an order-preserving FRDS in Cn> which means that 

77 < £ m Cd implies tp{t,u))r] < ip(t,uj)t; for allt>0 and u) G 0. 

Proof. This follows from Theorem [O] with G = G. □ 
Theorem 15.71 makes it possible to apply the general theory of monotone RDS 
(see [4] and also [21 [5]) to the class of sfde's considered. In particular it is possible 
to obtain the following results: 

• To provide transparent conditions which guarantee the existence of sto- 
chastic equilibria and a compact pull-back attractor (see, e.g., the general 
Theorem 3.5.1 in [4]). We recall (see [1]) that a random variable w : O j— > 
Cd is said to be an equilibrium (or fixed point, or stationary solution) of 
the RDS (1?, tp) if it is invariant under tp, i.e. if 

tp(t,uj)u(uj) = u(i?(t)w) a.s. for all t > 0. 

One can see that any equilibrium u(u) G Cd for ($, tp) has the form 
m(t,u;) = i>($(t)w), t G [— r, 0], where v(u>) is a random variable in D. 

• To describe the pull-back attractor for ($, tp) as a compact set lying be- 
tween two of its equilibria (see general Theorem 3.6.2 in [3]). 

• In the case when there exists a probability measure tt on the Borel a- 
algcbra of subsets Co, such that the law C(tp(t,oj)x) weakly converges to 
7r in Cd, the system ($, tp) has a random weak attractor A[uj) which is 
singleton, i.e. A(ui) = {i>(w)}, where the random variable v(uj) G Cd is 
an equilibrium (see the general Theorem 1 proved in [5]). We note that 
sufficient conditions for the existence of an invariant probability measure n 
and weak convergence of transition probabilities to ir for an sfde (monotone 
or not) have been established for example in [13] , [24] . |10] , and [TTj . 

Example 5.8 Consider the stochastic equations 

dx'(i) = (gHx 1 (t), . . . ,^(t))+ 5 U^ 1 (*-^i) I ■ • ■ ,x d {t-r d ))dt+m\x 1 )AW l (38) 

for i = l,...,d. We assume that g l , g\ and m l are smooth functions which 
arc globally Lipschitz. Under these conditions we can apply Proposition ^. 4f iii) 
to guarantee global well-posedness for (|3"8"]l . Moreover, one can see from The- 
orem 15.71 that equations (f38|) generate an order-preserving RDS in the space 
C = C([— r, 0]; R d ) with r = max^ ri provided that 

OXi 

and g\{x) is monotone, i.e. for every i = 1, . . . ,d we have that 

g\(x\...,x d )<gl(y\...,y d ) when ' //•*' . ./ 1 </. 
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The following example is a special case of Example 15.81 

Example 5.9 Let W be standard Brownian motion. Consider ID the retarded 
stochastic differential equation 

dx{t) = (f(x(t)) + g(x(t - l)))dt + a(x{t))dW(t), 

where /, g, a : R — > R arc Lipschitz, g is monotone and a is strictly positive. By 
Theorem 15.71 this equation generates an order-preserving RDS in C([— 1,0], R) 
Assume that the associated Markov semigroup on C([— 1,0], R) admits an in- 
variant (or stationary) measure (sufficient conditions are provided in |24] ) . In 
this case we can apply Theorem 1 [5] and conclude that the corresponding RDS 
has a unique equilibrium which is a weak random attractor. 

It is known (cf. [53]) that in the case a = the attractor of this system can 
contain multiple equilibria and also a periodic orbit. Thus we observe here that 
adding the noise term simplifies essentially long-time behavour of the system 
(see also [5] for some details). 

Example 5.10 (Stochastic biochemical control circuit) We consider the 
following system of Stratonovich stochastic equations 

dx 1 ^) = {g{L d x d t ) - a 1 x 1 (t))ch + a x ■ x x {t) o dW] , (39) 

dx°{t) = (Lj^x^ 1 - a 3 x ] {t))dt + o-j ■ x 3 (t) o dWf , j = 2,...,d. (40) 

Here as above "o" denotes Stratonovich integration, (jj are nonncgative and aj 
are positive constants, j = 1, . . . , d, and g : M + n> M + is a C 1 function such 
that 

< g(u) < au + b, and g'(u) > for every u > 

for some constants a and b. We also use the notation Xj t t{s) = Xj(t + s) for 
s G [— Tj, 0] and 

where fj,j : [— rj,0] t-> R is nondecreasing, ^j(-rj) = 0, (J,j(0) = 1, Mi( s ) > f° r 
s > —Tj. We denote r = maxj rj and equip (|39]l and (|40| with initial data 

x i (t) = C(t)>0, te[-r,0], i = l,...,d. (41) 

A deterministic version of this system was considered in [28], the stochastic 
non-retarded case was studied in [3] , see also [5] . 

Let £ e C([-r,0] : Rf ) by Proposition [23] a local solution 

x(t) = (x\t),...,x d (t)) 

exists on some interval [0, T(w)). By Theorem 13.41 fsee also Remark 13.21 and 
Corollarv l3.6[) we have that R^. is a forward invariant set, i.e. x(t) £ R^ for all 
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t G [0, r(w)). Applying Comparison Principle (see Theorem 14. 2 p we conclude 
that 

0<x(t)<x(t) for all t G [0,T(w)), (42) 
where x 4 (i) = (x 1 ^), . . . ,x d (t)) solves the following system of linear equations 

dx 1 (t) = (aL d x d ~ aix 1 (t) + b)dt + a t ■ x 1 (t) o dW? , (43) 

dx j (t) = (Lj-ixr 1 - UjX j (t))dt + aj ■ x 3 (t) o dWl , j = 2,...,d. (44) 

with initial data (|4"Tj) . The structure of (|43p and (j3U) allows us to solve these 
equations. Indeed, if we consider the drift part of the problem: 

dx j (t) = aj ■ x j (t) o dWj , x J '(0) = x, j = l,...,d., 

then ip 3 (t,x) = x exp{(TjW^ (t)} solves it. The spatial derivative of x) and 
its inverse are are independent of x and thus we can apply Proposition 12. 4f iii) 
to prove global existence of the solution x. Due to (|42|) this implies that the 
solution x(t) of (|31?|) . P0|) and (|4"Tj) docs not explode at finite time and thus 
equations fl39]) and (gUJ) generate an RDS (&,ip) in C+ = C([-r,0] : where 
r = max^ri. By Theorem 15.71 this RDS is order-preserving. By Comparison 
Theorem 14.21 this system is dominated from above by the affinc RDS (•&, <p a f) 
generated by (03) and (|44|) . 

Now we concentrate on the case a = (this means that g(u) is bounded). In 
this case we can construct an equilibrium v(to) = (u 1 (o;), . . . , v d (uj)) for (d, (f a f) 
by the formulas 

vi(u) = b / e ait - aiW ' dt, 

J — OO 

and 

w» = / Li-i^r 1 • e^'-^dt, j = 2, . . . ,d, 

J — oo 

where v 3 t {uj) :— v J ^(t + r)), r G [— r, 0]. Since (p(t,uj)x < (p a f(t,u>)x for every 
x G C+, it is easy to see that v{ui) is a super-equilibrium for ($, </?), i.e., 

< v(i3tij) a.s. for all t > 0. 

Thus by Theorem 3.5.1 g] the RDS (0,<p) has an equilibrium G R%. If 

g(0) > 0, this equilibrium is strongly positive. 

We can also show that in the case a = the RDS </?) possesses a random 
pull-back attractor in the universe of all tempered subsets of C([— r, 0] : R+). 
Indeed, due to the compactness property of the cocycle *& (see Theorem 15. 3j) it 
is sufficient to prove that possesses a bounded a absorbing set. This set 

can be constructed in the following way. 

Let v(lu) be the equilibrium for ($,ip a f) constructed above. One can see 
that in this case v\(u>) = Xv(u>) is a super-equilibrium for RDS ($, (/? a /) for 
every A > 1. One can also see that the top Lyapunov exponent for ($, ip a f) with 



24 



a = b = is negative. This implies v\(uj) is an absorbing super-equilibrium for 
(i?, (faf), i-e. for every tempered set D(uj) in C([— r, 0] : Rl) there is to{oj) that 

tp a f(t,'d-tOj)y('&-tU)) < v\(u), t>t D (u>), y e D. 

Since ( , d,ip a f) dominates tp), this implies that the interval 

[0, v x (cj)} = {ueC : 0<u< v x {uj)} 

is absorbing for </?). Therefore Theorem 1.8.1 [J] implies the existence of 
a pullback attractor which belongs to some interval of the form [ui(ui), u?,(ijj)\, 
where [iti(w) and v,2(u))] are two equilibria such that < u\(uS) < u%{u3) < v(lu). 
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